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Introduction 

In  this  paper  we  describe  and 
compare  methods  of  detecting 
short,  randomly  occurring  mes¬ 
sages.  The  emphasis  is  on  spread- 
spectrum  transmissions,  but  the 
energy  detection  results  also  apply 
to  simple  burst  communications, 
such  as  those  used  for  interrogation 
or  identification  among  network  ele¬ 
ments.1 

A  historical  account  of  spread- 
spectrum  systems  is  given  in 
Scholz.2  Some  early  papers  con¬ 
cerning  implementation  have  been 
reprinted3.  Dixon’s  book4  is  a  good 
introduction  to  spread-spectrum 
techniques.  Pickholtz  et  al.5  provide 
a  tutorial  treatment  of  the  theory.  A 
three-volume  set  by  Simon  et  al.6 
includes  a  chapter  on  the  detection 
of  low  probability  of  intercept  (LPI) 
signals.  Torrieri1  covers  secure 
communications  and  discusses  the 
intercept  problem.  Schleher7  in¬ 
cludes  a  discussion  of  the  detecta¬ 
bility  of  LPI  waveforms.  Wiley8  cov¬ 
ers  methods  of  intercepting  radar 
LPI  waveforms.  Dillard  and  Dillard9 
address  the  detection  of  spread- 
spectrum  signals.  Several  journal 
papers10-14  outline  techniques  for 
detecting  LPI  signals. 

Signal  Structures 

Figure  1  illustrates  the  use  of  no¬ 
tation  and  terminology.  The  duration 
and  bandwidth  of  the  transmission 
are  denoted  by  ^  and  Wt.  For 
hepped  time  or  frequency  transmis¬ 
sions,  the  duration  and  bandwidth  of 
the  frequency-hopped  pulses,  or 
data  symbols,  are  T2  and  W2.  In  a 
pure  frequency-hopped  (FH)  trans¬ 
mission,  the  data  symbol  is  a  simple 
pulse  (T2W2=1),  and  the  pulses  are 
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contiguous  in  time.  In  a  pure  time- 
hopped  (TH)  transmission,  the  data 
symbol  is  a  simple  pulse  and  all 
occupy  the  same  frequency  band 
(Wi  =  Wz).  The  duration  of  the  TH- 
interval  is  denoted  by  Tp.  Pseudo¬ 
noise  (PN)  is  a  wideband,  noise-like 
transmission.  A  pure  PN  transmis¬ 
sion  would  fill  the  area  T^Wi.  A 
hybrid  signal  is  a  combination  of 
two  or  more  spread-spectrum 
techniques.  The  most  common 
are  FH/PN,  FH/TH,  TH/PN,  and 
FH/TH/PN. 

The  number  of  data  symbols,  or 
pulses,  in  a  hopped  transmission 
is  denoted  by  b,  with  b  =  T-|/T2  in 
an  FH  or  an  FH/PN  transmission 
and  b  =  T,/Tp  in  a  TH,  TH/PN,  or 
FH/TH/PN  transmission.  The  sig¬ 
nal  energy  of  the  transmission  is  Ei 
and  the  energy  of  each  data  symbol 
is  E2  =  Ei/b. 

The  primary  detection  technique 
useful  against  weak  PN  and  burst 
signals  is  energy  detection.  If  the 
transmission  involves  hopped  sig¬ 
nals,  then  additional  detection  tech¬ 
niques  should  be  considered.  A 
wideband  energy  detector  that  inte¬ 
grates  energy  over  the  entire  trans¬ 
mission  is  often  the  best  choice. 

Energy  Detection 

Ideal  energy  detection  is  de¬ 
scribed  by  the  equation 

t 

V(t)=  |  y2(r)dr 

t-T  0) 

where  y(t)  is  the  input  voltage,  V(t) 
the  output  voltage,  and  T  the  inte¬ 
gration  time.  When  used  in  an  inte¬ 


grate-and-dump  (ID)  mode,  the  out¬ 
put  of  the  integrator  is  sampled  ev¬ 
ery  T  seconds  and  then  reset  to 
zero.  If  two  overlapping  integrators 
are  used  with  one  filter,  then  a  de¬ 
cision  is  made  every  T/2  seconds 
for  that  filter.  In  the  case  of  contin¬ 
uous  integration,  the  output  from  the 
integrator  at  any  time  t  is  the  energy 
integrated  over  the  time  interval 
(t-T,t).  Digital  integration  can  provide 
a  constant  weighting  with  time  when 
continuous  integration  is  desired, 
but  usually  with  a  small  degrada¬ 
tion.15 

Filter-Bank  Energy-Detection 
Systems 

A  filter-bank  energy-detection 
system  is  effective  when  the  fre¬ 
quency  band  of  each  transmission, 
such  as  a  burst,  varies  from  trans¬ 
mission  to  transmission.  It  is  also 
effective  against  hopped-frequency 
pulses.  Against  FH  and  FH  hybrid 
transmissions  where  the  possible 
frequencies  are  known,  a  bank  of 
bandpass  filters  can  be  used  to  cov¬ 
er  the  frequency  range,  matching 
the  filter  bandwidth  W  to  the  band¬ 
width  W2  of  the  data  symbols. 

In  the  simplest  filter-bank  system, 
each  filter  is  followed  by  an  energy 
detector  and  a  thresholder.  While 
the  system  is  modeled  this  way,  the 
implementation  may  be  more  com- 


Fig.  1  An  illustration  of  a  signal-parameter 
notation. 
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plicated.  For  example,  an  RF  spec¬ 
trum  analyzer  implemented  by  us¬ 
ing  an  acousto-optic  Bragg  cell  with 
an  exposed  film  strip  as  the  output 
can  function  as  a  continuously  in¬ 
tegrating  bank  of  energy  detectors, 
one  filter  for  each  of  a  number  of 
frequency  bands.16  An  example  of 
an  l+D  system  is  a  Bragg  cell  spec¬ 
trum  analyzer  that  uses  an  array  of 
photosensors  containing  photosites 
that  are  sampled  (read  out)  and  then 
cleared.17  Schleher7  describes  a 


Bragg  cell  system  that  uses  photo¬ 
diode  arrays  coupled  to  charge- 
coupled  devices  for  readouts.  Per¬ 
formance  standards  for  Bragg  cell 
receivers  are  introduced  in  Tsui  et 
al.18 

Energy  Detection  Calculations 
When  an  energy  detector  is  used 
to  integrate  the  energy  of  the  entire 
transmission,  the  integration  time 
should  be  T  =  1!  and  the  passband 
of  the  filter  should  match  the  trans- 
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mission  band,  with  bandwidth  W  = 
Wv  Calculations  for  pulse-detection 
systems  use  the  concept  of  dividing 
the  time-frequency  space  under 
surveillance  into  cells  of  energy  in¬ 
tegration.  This  defines  a  grid  of  cells, 
each  cell  having  area  TW  *  T2W2. 
An  ideal  fit  of  the  grid  to  the  signal 
results  in  only  two  kinds  of  integra¬ 
tion  cells,  cells  containing  signal 
energy  and  cells  empty  except  for 
noise. 

To  analyze  energy  detection,  we 
assume  that  the  noise  at  the  input  to 
the  energy  detector  is  a  zero-mean, 
stationary,  Gaussian  random  pro¬ 
cess  that  has  a  flat,  band-limited 
power  spectral  density.  For  large 
TW,  we  assume  that  an  ideal  rectan- 
gular-passband  filter  of  bandwidth 
W  precedes  the  energy  detector. 
When  TW  is  in  the  region  of  unity, 
the  received  signal  energy  depends 
on  the  pulse  shape  and  filter  shape. 
The  normalized  decision  statistic  V 
is  related  to  the  integrated  energy  V 
by  V’  =  2V/N01,  where  N01  denotes 
the  one-sided  noise  power  spectral 
density. 

The  thresholder  compares  the 
value  of  V’  with  a  threshold  K,  and 
the  decision  is  "signal  present”  if  Y 
>  K,  and  “noise  only”  if  V  <  K.  For 
continuous  integration,  V’  is  the 
value  of  the  normalized  integrated 
energy 

t 

v'(t)  ‘  KvT  |  yI(,)dr 

t-T  (2) 

at  any  particular  instant  t  For  l+D 
detection,  the  decision  times  are  at 
instants  T  seconds  apart  The  false- 
alarm  probability  per  decision  is 

Qf  =  Prob[  V'>K  at  time  t]  no  signal 
energy  in  (t-T,  t)]  (3) 


and  the  detection  probability  per 
decision  is 

Qo  =  Prob[  V’  >  K  at  time  tj  signal 
energy  Es  in  (t-T,  t)].  (4) 

The  vertical  bar  |  indicates  that  the 
probability  is  conditional  on  the 

(Continued  on  page  140) 
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event  following  the  bar. 

The  distribution  function  of  the 
test  statistic  V  is  closely  approxi¬ 
mated  by  the  chi-square  distribution 
with  2TW  degrees  of  freedom  in  the 
noise  case  and  by  the  noncentral 
chi-square  distribution  with  2TW 
degrees  of  freedom  and  noncen¬ 
trality  parameter  A.  =  2Es/Noi  in  the 
signal  case.1-10  For  TW  =  1,  the  dis¬ 
tribution  is  equivalent  to  the  square- 
law  form  of  the  Rice  distribution, 
which  describes  envelope  detec¬ 


tion  of  a  single  pulse.  When  TW  is 
any  integer,  the  equation  for  Qo  is 
the  generalized  Q-function,  Qo  = 
QTW( V2S.-/K),  where  S  = 
Es/N0i  is  the  signal-to-noise  ratio 
(SNR).10 

The  value  of  the  bias  level,  the 
normalized  threshold  Ko  =  K/2,  that 
provides  the  desired  Qf  is  found  in 
Pachares’  table19  for  TW  = 

1.2..  ...150,  and  for  QF  =  10‘p,  P  = 

1.2..  ...12.  A  technique  outlined  by 
Helstrom20  can  be  used  to  approx¬ 


imate  Ko  for  moderately  large  values 
of  TW,  such  as,  for  TW  >  20. 

Robertson21,  Dillard22,  and  Pari23 
give  recursive  formulas  for  comput¬ 
ing  Qo.  Helstrom20  gives  techniques 
for  an  approximate  evaluation.  This 
approximation  and  others  are  dis¬ 
cussed  below. 

Energy  Detection  Approximations 

For  large  values  of  TW,  the  central 
limit  theorem  can  be  used  to  obtain 
normal  approximations  to  QF  and 
Qo.10  In  terms  of  the  normal  cumu¬ 
lative  distribution  function 

y 

F(y)  =  vkr  f  e"2/2dz 

(5) 

the  approximations  are 

Qf~1  -F[(K0-TW)/v/TW]  (6) 


and 

QD  =1  -  F  [(Ko-TW-S)/  \/TW+2Sj. 

(7) 


By  solving  Equation  6  for  the  nor¬ 
malized  threshold  Ko  we  obtain 

Ko  *»  F'^I-QfJn/tW  +  TW,  (8) 

where  F'^x)  is  equal  to  the  variate 
y  such  that  F(y)  =  x.  Substitution  of 
Equation  8  into  Equation  7  yields 

Qo  =1  -  FffF-^l-Qp)  -  S/VlW 
Vi  +  2S/TW).  (9) 


The  approximation  in  Equation  9  is 
called  the  "full  normal  approxima¬ 
tion”.  A  further  approximation  valid 
for  S«  2TW  is 

Q0  =  1  -  F[F-1(1  -QF)-S/VTWj. 

(10) 


The  approximations  given  by 
Equation  9  and  Equation  10  for 
energy  detection  require  the  prod¬ 
uct  TW  to  be  large,  but  their  accu¬ 
racy  also  depends  on  the  values  of 

(Continued  on  page  144] 
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HEMTs 
1  TO  30  GHZ 
ULTRA  LOW  NOISE 

Fujitsu  now  offers  a  complete 
family  of  ultra  low  noise  quarter 
micron  HEMT  devices  for  use 
from  L  through  Ka  band.  These 
devices  are  both  rugged  and  reli¬ 
able.  Space  qualified  units  can  be 
provided.  Uses  include  receiver 
front  ends  for  communication, 
radar,  and  satellite  systems. 

Typical  Performance 


Device 

Test  Freq. 
(GHz) 

N.F. 

(dB) 

Gas 

(dB) 

Gate  Width 
(Microns) 

FHC30 

4 

0.35 

14.5 

280 

FHC31 

4 

0.50 

14.5 

280 

FHX04 

12 

0.75 

10.5 

200 

FHX05 

12 

0.90 

10.5 

200 

FHX06 

12 

_ 

12.0 

200 

FHR02 

18 

1.0 

9.0 

200 

FHR10 

18 

1.0 

9.5 

100 

For  full  details  contact  Fujitsu  Microelectronics,  Inc. 
Microwave  and  Optoelectronics  Division  ‘ 

3330  Scott  Blvd.,  Santa  Clara,  CA  95054-3101,  (408)  562-1550 
TELEX:  278868  FMMO  UR,  FAX:  (408)  727-3194 

fuJitsu 

Th*  glflbj]  computer  &.  communications  company. 
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Qf  and  Q&  Equation  7  provides  a 
good  approximation  to  Q0  when  the 
value  used  for  Ko  is  accurate,  even 
for  relatively  small  values  of  TW, 
such  as,  TW=20  or  larger.  Equations 
9  and  10  are  in  error  primarily  be¬ 
cause  the  approximation  to  Ko  in 
Equation  8  is  in  error. 

If  Ko  is  known,  we  can  use  Equa¬ 
tion  7  to  estimate  the  SNR  required 
to  achieve  a  given  value  of  Qo.  The 
estimate  is 

S  =  K  o  -  TW+-/2  -  y\/2Ka-TW+72 

(11) 


where  y  =  F'T(1-QD).9  When  Qq  = 
0.5,  then  7  =  0,  and  the  estimate  is 

S(0.5)  =  Ko  -  TW.  (12) 


Equation  11  is  obtained  by  solving 
the  quadratic  equation  K0-TW- 
S=7\/TW+2S  for  S.  Note  from 
Equation  12  that  the  difference  Ko  - 
TW  is  an  estimate  of  the  required 
SNR  for  Qd=0.5. 

The  two  normal  approximation 
Equations  9  and  1 0  and  the  “hybrid" 
approximation  using  Equation  7  and 
the  true  Ko  are  compared  with  the 
“exact”  results  (using  the  noncen¬ 
tral  chi-square  distribution)  in  Fig¬ 
ure  2  for  Qf=10~12.  The  hybrid  ap¬ 
proximation  gives  results  equal  to 
the  exact  results  except  for  small 


TW,  where  none  of  the  approxima¬ 
tions  are  reasonable.  Therefore,  if 
the  true  values  of  Kq  are  available, 
there  is  no  need  to  use  Equation  9 
or  10.  Evaluation  of  Equation  7  using 
Hasting’s  approximations24  is  much 
simpler  than  evaluation  of  the  non- 
central  chi-square  distribution  us¬ 
ing,21*23  rather  than  the  approxima¬ 
tion  using  Helstrom's  technique.20 

False-Alarm  Rate 

The  false-alarm  rate  (FAR)  for  a 
single-filter  energy-detection  sys¬ 
tem  with  simple  thresholding  is  the 
product  of  the  decision  rate  and  the 
false-alarm  probability  Qf.  For  l+D 
detection,  the  decision  rate  is  1/T; 
for  overlapping  l+D  detection,  it  is 
2/1.  The  FAR  for  the  filter-bank  sys¬ 
tem  with  simple  thresholding  is 

FAR  =  (number  of  filters)  x  (decision 
rate  per  filter)  x  Qp 

For  l+D  detection,  we  then  have 
FAR  =  NQp/T,  where  N  is  the  number 
of  filters.  In  Figure  3,  the  FAR  is  held 
constant  with  changing  TW  by  fixing 
the  value  of  Qp/TW.  Note  the  mis¬ 
leading  fact  that  Qo  improves  as  TW 
increases  when  computations  stray 
into  the  area  where  Q0  is  less  than 
10  Qf-  When  integration  is  over  the 
entire  transmission,  the  message 
detection  probability  is  Pd  =  Qd- 
When  integration  is  over  individual 
data  symbols,  the  message  detec¬ 
tion  probability  is  a  function  of  Qd 


and  b.  Figure  3  applies  to  both 
cases.  Curves  for  other  false-alarm 
rates  are  given  in  Reference  9. 

Splitting  Loss 

For  l+D  detection,  matching  the 
integration  time  to  the  signal  dura¬ 
tion  is  optimum  only  if  the  integra¬ 
tion  interval  coincides  with  the  sig¬ 
nal.  In  practice,  the  signal  is  likely  to 
be  split  into  two  intervals  of  integra¬ 
tion.  Figure  4  gives  an  example  of 
the  effect  on  detection  probability, 
for  TW  ~  1.  The  average  detection 
probability  is  1  -  [1  -QD(xE)][1  -QD(E 
-  x  E)]  averaged  over  x,  where  x  is 
uniformly  distributed  between  0  and 
0.5.  Figure  5  gives  the  loss  Lt  corres¬ 
ponding  to  the  average  detection 
probability  for  several  combinations 
of  TW  and  Qp 

Mismatched  Integration  Time 

If  the  integration  interval  signifi¬ 
cantly  exceeds  the  duration  of  the 
signal,  a  loss  cesults  from  the  un¬ 
necessary  integration  of  additional 
noise  power.  A  small  offset  to  this 
loss  is  the  decreased  splitting  loss. 
If  the  integration  time  is  significantly 
less  than  the  signal  duration,  then 
detection  involves  “repeated  obser¬ 
vations".  For  message  duration  T,, 
one  computes  QD  for  Es  =  EiT/T, 
and  lets  the  message  detection 
probability  equal 

Pd*1 -(1 -Qo)Tl/T  (13) 

[Continued  on  page  146] 


Fig.  2  Comparison  of  approximations  with  exact  distribution. 


Fig.  3  Detection  probability  per  integration  for  a  fixed  FAR. 
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Energy  Integration 
over  Each  Data  Symbol 

Consider  a  hopped  transmission 
consisting  of  b  data  symbols.  The 
data  symbols  can  be  simple  pulses 
(T2W2  =  1;  as  in  FH,  TH,  or  FH/TH), 
or  can  have  a  PN  microstructure 
(T2W2»1;  as  in  FH/PN,  TH/PN,  or 
FH/TH/PN).  The  detection  system 
is  modeled  by  a  bank  of  N  fiiter- 
detector-thresholder  units,  where  N 
=  1  for  TH  or  TH/PN.  By  using  the 
grid  concept,  there  are  b  cells  con¬ 
taining  signal  energy  and  m  noise- 
only  cells  in  the  period  of  length  T1t 
where 


m  =  N 


T2 


b 


(14) 


If  a  signal  decision  in  one  of  the 
m  empty  cells  leads  to  an  examina¬ 
tion  of  the  message  in  a  recognition 
stage,  then  Pd  is  the  probability  that 
a  threshold-crossing  occurs  during 
the  message  arrival  time.  For  an 
ideal  grid  fit,  we  have 

Pd  =  1  -  (1  -  Qo)b  (1  -  Or)" 

(15) 


Fig.  4  Example  of  splitting  signal  energy  into  two  time  intervals  of  energy  integration. 


Splitting  loss  applies  only  to  the  first 
and  last  integration  of  signal  energy. 

Figure  6  gives  an  example  of  the 
loss  for  both  kinds  of  mismatching, 
for  a  fixed  message  size  (TV W,  =  1 04) 


and  for  varying  T,  with  W  =  Wv  Re¬ 
call  that  Qf  varies  with  T  for  a  fixed 
FAR.  Figure  7  shows  the  average 
loss,  which  varies  only  slightly  with 
TW,  Qf,  and  QD. 


where  Qd  is  the  probability  of  de¬ 
tecting  a  data  symbol.  The  approx¬ 
imation 

Pd~1-U-Qo)b  06) 

[Continued  on  page  148[ 


Fig.  5  Loss  as  a  function  of  the  probability  that  the  signal  is 
detected  in  at  least  one  cell  averaged  over  all  splitting 
proportions. 


Fig.  6  Illustration  of  result  of  mismatching  the  integration  interval 
T  with  the  signal  duration  Tv 


146 


MICROWAVE  JOURNAL  •  SEPTEMBER  1989 


[From  page  146]  DILLARD 


is  usually  valid.  Note,  however,  that 
if  Qq  is  close  to  Qf,  (1 5)  gives  Pd  = 
bmQp  while  Equation  16  gives  Pd  *= 
bQh  and  these  can  differ  by  orders 
of  magnitude. 

The  results  above  assume  that 
each  data  symbol  is  contained  en¬ 
tirely  in  some  time  bandwidth  cell  of 
integration.  The  loss  L(  resulting 
from  a  random  splitting  of  the  signal 
into  two  adjacent  passbands  can  be 
taken  into  account  when  estimating 
the  received  energy  E2.  A  loss  Lt 
occurs  for  l+D  detection  when  the 
signal  is  split  randomly  into  two  in¬ 
tegration  intervals.  For  PN  data  sym¬ 
bols,  losses  Lf  and  L,  are  approxi¬ 
mately  equal  to  each  other  and  are 
additive.9 

Binary  Moving-Window  Detection 

A  binary  moving-window  detec¬ 
tor  (BMWD)25  implements  the  re¬ 
quirement  that  at  least  k  out  of  b 
pulses  be  detected  for  a  “signal  pre¬ 
sent”  decision  to  be  made.  In  the 
applications  here,  some  of  the  de¬ 
tections  result  from  noise  in  other¬ 
wise  empty  time-bandwidth  cells. 


First  consider  the  use  of  a  BMWD 
system  against  TH  or  TH/PN  sig¬ 
nals.  An  energy  detector  (with  inte¬ 
gration  time  T  =  T2)  followed  by  a 
thresholder  provides  the  input  to  the 
BMWD.  Every  T2  seconds  the  input 
is  a  0  or  a  1.  (A  run-suppressor  unit 
is  also  needed,  because  the  inte¬ 
gration  interval  is  not  synchronized 
with  pulse  times.)  The  binary  inte¬ 
gration  time  is  TM,  usually  the  mes¬ 
sage  duration,  T,.  At  each  time  jT2, 
the  BMWD  computes  the  test  statis¬ 
tic 

i 

si-  X  y* 

k  =  j-M+1  (17) 

and  compares  it  with  a  threshold  Km- 
The  window  length  is  M,  where  M 
=  Tm/T2.  An  alarm  occurs  when  the 
test  statistic  Sj  reaches  the  threshold 
Km  in  an  upward  direction. 

If  synchronization  with  the  hop¬ 
ping  intervals  of  a  TH  or  TH/PN 
transmission  were  possible,  the 


TEM  filter  with  a  response  identical  to  that  of  a 
“parallel-coupled  stripline"?  And  that  an  INTERDIGITAL 
filter’s  high-pass  design  in  Richard's  tangent  variable 
results  in  symmetrical  skirts  and  odd-harmonic 
passbands  (often  “cleaned  up”  by  cascading  a  low-pass 
filter)?  Delta’s  INTERDIGITAL  design  is  very  predictable, 
and  the  units  are  easily  fabricated  and  tuned. 

Unique  filter  requirements?  Talk  to  Delta! 


840  Via  Alondra 
Camarillo,  CA  93010 

(805)  987-6892 
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knowledge  that  one  and  only  one 
data  symbol  occurs  in  each  TH  in¬ 
terval  could  be  used  by  employing 
a  binary-tapped  delay-line  and  an 
OR-gate.  The  output  of  the  OR-gate 
every  Tp  seconds  would  then  be  y* 
=  0  if  every  slot  in  the  TH-interval 
yields  0,  and  would  be  1  otherwise. 
Computations  for  this  synchronized 
case  give  a  close  upper  boundary 
to  the  detection  capability  of  pulse- 
detection  systems  against  a  TH  or 
TH/PN  transmission,  and  are  much 
simpler  than  for  the  unsynchronized 
case. 

Against  FH,  FH/PN,  FH/TH,  and 
FH/TH/PN  transmissions,  a  filter- 
bank  BMWD  system  can  use  the 
knowledge  that  only  one  frequency 
slot  is  occupied  at  any  instant  by 
employing  an  OR-gate  to  allow,  at 
most  a  single  “1”  to  enter  the  win¬ 
dow  every  T  seconds.  The  synchro¬ 
nization  assumption  also  can  be  ap¬ 
plied  to  filter-bank  BMWD  systems 
for  detecting  FH/TH  and  FH/TH/ 
PN  transmissions.  In  the  latter  case, 
the  BMWD  accepts  a  single  1  or  0 
from  each  block  of  NTP/T2  cells. 

The  probability,  at  any  time  jT2  that 
the  binary  number  1  enters  the 
BMWD  is  denoted  by  p0  when  only 
noise  is  present  during  the  previous 
T2  seconds,  and  is  denoted  by  p-i 
when  a  signal  is  present  (in  one  of 
the  N  passbands).  For  the  single¬ 
filter  system  (unsynchronized  with 
the  TH-interval),  the  formulas  for  p0 
and  p,  are  simply  p0  =  Qf  and  pi  = 
Qo-  For  the  filter-bank  system  with 
an  OR-gate  applied  to  the  N  filter- 
detector-thresholder  inputs,  the  for¬ 
mulas  are 


[Continued  on  page  ISO] 


Fig.  7  Loss  resulting  from  mismatching  the 
integration  interval  with  the  signal  duration 
averaged  over  a  range  ot  parameter 
combinations. 
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Po-1  -(1'Qf)n  (18) 

and 

p,  =  1  -(I-Q^'HI'Qd)  (19) 
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At  each  time  jT 2,  the  BMWD  com¬ 
putes  Sj,  given  by  Equation  17,  and 
compares  it  with  the  threshold  Km, 
where  M  =  T1/T2  is  the  window 
length.  The  system  alarms  when  the 
sum  Sj  reaches  the  threshold  Km  in 
an  upward  direction,  such  as,  when 
Sj-i  =  Km-1  and  S,  =  KM.  The  prob¬ 
ability  of  this  joint  event  when  no 
signal  is  present  is 

/  M-1  i 

Pt  -  (  K  ]  PoKM(1  ~  po)u'Kw*’ 

(20) 

The  FAR,  in  alarms  per  second,  is 
FAR  =  PWT2. 

The  probability  of  a  detection  by 
the  BMWD  is  approximately  equal 
to  (and  bounded  below  by)  the  prob¬ 
ability  that  there  are  at  least  KM  Is 
in  the  window  when  the  window  is 
aligned  in  time  with  the  transmis¬ 
sion.  This  probability  is11 


M 

min(i.b) 

X 

i~kM  j 

=  max  (0,i-M+b) 

/  M  -b  \ 

(i  ) 

1  i-j )  p’,p”" 

•  (1  -  p,)b-i(1  -  po)M_id-b.  (21) 


For  FH  or  FH/PN,  with  M  =  b,  Equa¬ 
tion  21  becomes 

b 

Pd  =  S  (  !  )  p1‘n  ’  Pl)^ 

i  =  Km  (22) 

Equation  22  also  applies  to  systems 
assumed  to  have  an  OR  function 
synchronized  with  theTH-interval,  if 
we  use 

p-  =  1  -(1  -  Qf)ntp/T2'1(1-Qo).  (23) 

The  window  length  is  M=b,  and  the 
FAR  becomes  FAR=P-[7Tp,  where 
the  computation  of  PT  uses 

[Continued  on  page  152] 


Fig.  8  A  comparison  of  pulse  detection  systems  for  TH/PN  transmissions. 
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Po  =  1  -  (1  -  Qf)NTp/T2‘  (24) 


For  the  unsynchronized  systems 
used  against  signals  involving  TH, 
the  use  of  Equation  22  with  Equation 
23  serves  as  a  good  approximation 
to  Equation  21 .  For  fixed  FAR  and  E2, 
it  provides  an  upper  boundary  to 
Equation  21  when  Equation  24  is 
used  with  M  =  b  and  FAR  =  Pt/Tp, 
and  a  lower  boundary  to  Equation 
21  when  Equation  18  is  used  with  M 
=  bTp/T2  and  FAR  =  Pt/T2.  The 
value  of  Qr  and  therefore  of  pi,  is 
different  for  the  two  calculation 
methods. 

Comparisons 
of  Detection  Methods 

Against  hopped  transmissions, 
the  dedicated  interceptor  must  de¬ 
termine  whether  integration  over  the 
transmission  (IOT)  or  integration 
over  individual  data  symbols  (pulse 
detection)  is  better.  If  the  pulse-de¬ 
tection  approach  is  competitive, 
then  additional  signal-processing 
techniques,  such  as  BMWD  should 


be  considered.  Generally,  IOT  is  bet¬ 
ter  when  signal  energy  densely  oc¬ 
cupies  the  transmission’s  time- 
bandwidth  area  and  pulse-detec¬ 
tion  is  better  when  pulses  are 
sparsely  distributed. 

Against  TH/PN  signals,  the  size 
T2W2  of  the  data  symbols  strongly 
affects  pulse-detection  perform¬ 
ance,  as  illustrated  in  Figure  8.  Note 
that  IOT  becomes  better  than  pulse 
detection  as  T2W2  increases. 

Figure  9  compares  systems  for 
two  FH  cases,  one  where  only  ten 
frequencies  are  used  and  one 
where  104  are  used.  Note  that  for 
IQ4  frequencies,  IOT  becomes  bet¬ 
ter  than  simple  pulse  detection 
when  b  exceeds  3000,  and  better 
than  BMWD  with  OR-gating  when  b 
exceeds  104. 

The  results  shown  in  Figure  8  for 
TH/PN  and  in  Figure  10  for  FH  are 
reflected  in  Figure  10,  which  con¬ 
siders  IOT  detection  vs.  pulse  detec¬ 
tion  with  BMWD.  The  upper  left  area 
corresponds  to  sparse  signal  occu¬ 
pancy,  and  the  area  under  the 
curves  corresponds  to  dense  occu- 

[Continued  on  page  1 54] 


Fig,  9  A  comparison  of  detection  methods  for  pure  FH,  N  filters. 
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Fig.  1 0  For  combinations  of  T,  W,  and  b  that  fall  above  the  curve  corresponding  to  the 
selected  pulse  size,  pulse  detection  using  a  BMWD  is  the  better  detection  method. 
Below  the  curve,  integration  over  the  transmission  is  better. 


pancy.  For  intermediate  values  of 
signal  density,  the  designer  of  the 
detection  system  should  use  the  for¬ 
mulas  given  here  or  use  graphs9  to 
determine  whether  to  integrate  the 
signal  energy  over  the  entire  mes¬ 
sage  or  to  use  a  pulse-detection 
system.  ■ 
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